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Nuclear double B~ -decay with two neutrinos is an important decay mode for some unstable nuclei. Based 
on the available experimental data of nuclear double B~ -decay, we propose that there is a law between the 
logarithm of double B~ -decay half-lives and the reciprocal of the decay energy. The physics behind the law is 
discussed and it is found that this is associated with the universal properties of the weak interaction. This double 
B~ -decay law is similar to the famous Geiger-Nuttall law of «-decay. The law is applied to predictions of the 
nuclear double B~ -decay half-lives for six even-even nuclei from Z = 84 to Z = 98 and we found that *’Th 
is very interesting for future experiments. The branching ratios between double B~ -decay and «-decay are also 
estimated for the six even-even nuclei and this is useful for future experimental search of new emitters of double 


B~ -decay. 
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I. INTRODUCTION 


Researches on the B decay of nuclei lead to some important 
developments in both nuclear physics and particle physics. In 
order to explain the missing energy in B-decay, Pauli pro- 
posed in 1930s that there is a new particle, neutrino, who has 
carried away the missing energy. Fermi founded the basic the- 
ory of nuclear -decay where the neutrino is involved in the 
theory although it is not discovered in 1930s. In 1950s, Lee 
and Yang proposed that parity can not conserve for a weak 
process such as B-decay [1]. Wu and her collaborators carried 
out the B-decay experiment with polarized Co nuclei and 
confirmed that parity symmetry was broken in the B- decay 
process [2]. This brought to important impact to the devel- 
opment of physics. For double B-decay, it was first predicted 
by Goeppert-Mayer in 1935 [3] and it was later observed by 
the experiments of geochemistry, radiochemistry and nuclear 
physics [4-15]. Wu also gave good advice on the direct obser- 
vation of rare double -decay and finally the double B-decay 
of *?Se was directly observed by Elliott, Hahn, and Moe [4]. 
To date, eleven nuclei are observed to have double B~ -decay 
[10-13]. Although the data of double B-decay was accu- 
mulated, no simple law among the data was found. This is 
very different from the researches of other radioactivity such 
as «-decay and cluster-radioactivity where a unified formula 
between the half-lives and decay energies is proposed [16]. 
A new Geiger-Nuttall law [17] has been also proposed very 
recently and the influence of quantum numbers on «-decay 
half-lives is included in the formula [17, 18]. For numerical 
calculations of double -decay half-lives, much of the theo- 
retical research contains complicated calculations with long 
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codes in computers and can only be done by professional the- 
oretical physicists. Of course some modern shell model cal- 
culations are still not available for heavy nuclei, such as 238U, 
due to the limit of the computational ability of computers. 
Due to the computational complexity of double B-decay half- 
lives, double -decay processes have not often been discussed 
in standard textbooks for undergraduate students and graduate 
students, although double 6-decay has been clearly observed 
for many years. This is very different from both a-decay and 
complex cluster radioactivity, such as '4C, where their half- 
lives can be easily calculated by simple formulas with three 
or four parameters. Therefore a simple and accurate way of 
calculating of double B-decay half-lives is needed. A simple 
and accurate way, and it is also useful for experimental physi- 
cists to estimate the decay half-live before an experiment of 
double (-decay and to analyze the data after the experiment. 


II. NUMERICAL RESULTS AND DISCUSSION 


Recently we make a systematic analysis on double B~- 
decay data and propose a new systematic law between the 
decay half-lives and the decay energy [19]. The law is sim- 
ple and accurate for the ground state transition of double B~- 
decay between parent nuclei and daughter nuclei [19]. The 
law, without any extra adjustments, also works well for the 
double ~-decay between the ground state of parent nuclei 
and the first excited state 0* state of daughter nuclei [19]. 
Because the systematic law is proposed by analyzing the ex- 
perimental data of the ground-state transitions of eleven even- 
even nuclei from “8Ca to PSU for double 8 ~-decay, the effect 
of the weak interaction, the correction of the Coulomb inter- 
action and the leading term of nuclear structure are taken into 
account in the law. The systematic law of double B-decay 
half-lives with two neutrinos is as follows [19] 


lg Ty/2(Ey) = (a — 2lg(2m 2/137) + 5)/ Q28; Cd) 
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where T}/2 is the double B~-decay half-life and its unit is 
Ey (10!8 years), Q2g is the double B~-decay energy (MeV), 
Z is the charge number of the parent nucleus [19]. The first 
term of this equation corresponds to the impact of the weak 
interaction and the constant a equals to 5.843 by fitting the ex- 
perimental double B~ -decay data of ground state transitions 
for eleven even-even nuclei. The second term of equation (1) 
corresponds to the effect of the Coulomb potential for double 
fs -decay. The third term, however, is related to the shell ef- 
fect and S' is chosen to be S = 2 when the neutron number 
of the parent nucleus is a magic number, S = 0 when the 
neutron number of the parent nucleus is not a magic number. 
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Fig. 1. (Color online)Variation of the logarithms of calculated and 
experimental double -decay half-lives with the decay energy for 
eleven even-even nuclei from “Ca (Z = 20) to SU (Z = 92). It 
is interesting to note that the half-lives also depend on the charge 
number from Eq. (1). It is seen that the half-lives are very sensitive 
to the decay energies. 


The variation of the logarithms of double B-decay half- 
lives with the decay energy is drawn in Fig. | for eleven even- 
even nuclei from “8Ca (Z = 20) to PSU (Z = 92). In Fig. 1, 
the X-axis is the decay energy of ground-state transitions and 
the Y-axis is the logarithm of double f-decay half-lives. Al- 
though the decay energy vary in a narrow range with several 
MeV, the half-lives varies in a very wide range from Ey (1018 
years) to 10° Ey. This demonstrates the important effect of 
decay energies on half-lives. The higher decay energy usu- 
ally leads to a shorter half-life for the nuclei on an isotopic 
chain. The calculated results with Eq. (1) reasonably agree 
with the available data and the number of adjusting parame- 
ters in Eq. (1) is the minimum, as compared to other theoret- 
ical calculations. The average deviation of half-lives of the 
eleven nuclei is a factor of three, which is good for a very 
simple formula of Eq. (1). 

In the following context, we predict the double -decay 
half-lives of six radioactive nuclei from 7!°Po (Z = 84) to 
256Cf (Z = 98) by using Eq. (1). The calculated results of the 
six nuclei are listed in Table 1 and drawn in Fig. 2. For the 
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Fig. 2. Variation of the logarithms of calculated double f-decay 
half-lives with the decay energy for six even-even nuclei from 7'°Po 
(Z = 84) to °C£(Z = 98). It is interesting to note that the half- 
lives are also dependent on the charge numbers from Eq. (1). It is 
seen that the half-lives are very sensitive to the decay energies. The 
double B-decay half-life of *’Th (Z = 90) is approximately 3 x 10° 
Ey and this shows that “Th is a very good candidate for experimen- 
tal observation of double ®-decay because it is rich in nature and it 
has a very long «-decay half-life (Tho = 1.4 x 10*° years). 


six nuclei, none of their neutron numbers are magic numbers, 
which means S$ = 0. 


In Table 1, the first column denotes the parent nucleus 
and the second column represents the experimental dou- 
ble B-decay energy for the ground-state transition between 
the parent nucleus and the daughter nucleus where Q2g = 
M(A, Z) — M(A,Z — 2) [12, 13]. The calculated double 
B-decay half-lives (T; /2(theor.) ) are presented in the third 
column of the table and the units of the double f-decay half- 
lives are Ey (1018 years). We also listed the experimental 
a-decay half-lives (Ti/9(expt.)) of the nuclei in column 4 
for comparison with the data from references [12, 13]. The 
fifth column represents the branching ratio between double 
B~ -decay and -decay (T{/,(expt.)/ Tı /2(theor.)). 


It is seen from Table | that the shortest double B-decay 
half-life is 8.36 x 10? Ey for *°Cf and the longest is 3.64 x 
1013 Ey for 22°Rn. They correspond to the highest double B- 
decay energy and the lowest double B-decay energy, respec- 
tively. This shows that the decay half-life strongly depends 
on the decay energy. Although 7°°Cf and ?!°Po have shorter 
half-lives for double 6-decay, they are not suitable for future 
double -decay experiments because their «-decay half-lives 
are too short. ??°Rn and 7°4Cf are also not suitable for fu- 
ture search of new emitters of double B-decay because their 
double $-decay half-lives are too long and their «-decay half- 
lives are also too short. Only ??Ra and ??Th are interest- 
ing for future double -decay experiments due to their longer 
x-decay half-lives. Especially ”?Th which has a very long 
a-decay half-life (1.4 x 101° years) and shorter double ß- 


$20510-2 


HALF-LIVES OF DOUBLE B- -DECAY WITH... 


Nucl. Sci. Tech. 26, $20510 (2015) 


TABLE 1. Theoretical double B~ -decay half-lives (T,/2(theor.)) for some heavy double $~ -decay candidates where the experimental «- 
decay half-lives (T2 (expt.)) are listed for comparison. The double ßB-decay energy and the branching ratio between the double B~ -decay 


and the «-decay are also compared. The units of the double 8-decay half-lives are Ey (101° years). The experimental «-decay half-lives and 
double B-decay energies (MeV) in the table are from Refs. [12] and [13] respectively 


Nuclei Q2p(MeV) Ty /2(theor.)(Ey) 
216 Do 1.531 1.13 x 103 
220Rn 0.343 3.64 x 101° 
226Ra 0.473 6.19 x 10° 
ai 0.838 3.19 x 10° 
2540F 0.438 2.30 x 101° 
250 CF 1.5534 ” 8.36 x 10° 


a 232 


Ta (expt) 


145 (ms) 4.07 x 107° 
55.6 (s) 4.84 x 10738 
1600 (y) 2.58 x 10775 
1.4 x 10'° (y) 4.39 x 107+ 
60.5 (d) 7.21 x 107%° 
12.3 (m) 2.80 x 10776 


Th will be very interesting for future experimental observation of double ®-decay due to both its shorter double f-decay half-life and the longest 


o-decay half-life in six nuclei. ?3?Th is also very rich in nature and this makes the experiment easier and cheaper as compared to other isotopes. 
b For 256Cf, the estimated mass is used for the calculation of the double B-decay energies (MeV) [12] and it is denoted with a # symbol. 


decay half-life (3.19 x 10° Ey), this leads to its largest branch- 
ing ratio between double -decay and a-decay. This clearly 
shows that ?3?Th is very interesting for future double ®-decay 
experiments. 7*?Th is also rich in nature and this makes the 
experiments easier and cheaper as compared to other isotopes. 

The variation of the logarithms of double B-decay half- 
lives with the decay energy is also drawn in Fig. 2 for six 
even-even nuclei from 7!°Po to 7°°Cf. We can see again 
that the logarithms of double -decay half-lives are inversely 
proportional to the decay energies. Usually a higher decay 
energy corresponds to a shorter double (-decay half-life and 
of course there is also a dependence on the charge number of 
parent nuclei. 

Finally it is interesting to make a brief discussion about 
the law of double B-decay half-lives (equation (1)). For the 
systematic law of double $-decay half-lives, the first term in 
the law depends on the decay energy because the weak inter- 
action is universal for natural decay processes and the total 
effect from the weak interaction is not very sensitive to the 
change of nucleon numbers such as charge number. This is 
different from that of -decay where the total effect of the 
repulsive Coulomb potential is directly related to the charge 
number of the nuclei and «-decay occurs for ground states of 
medium and heavy nuclei [17, 18]. 
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